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We examine the nature of topological currents in black phosphorus when its inversion symmetry 
is deliberately broken. Here, the conduction and valence band edges are located at the F point of 
the rectangular Brillouin zone, and they exhibit strong anisotropy along its two crystal axes. We 
will show below that these salient features lead to a linear transverse neutral topological currents, 
accompanied also by a non-linear transverse charge current at the Fermi surface. These topological 
currents are maximal when the in-plane electric field is applied along the zigzag crystal axes, but 
zero along the armchair direction. 

PACS numbers: 72.80.Vp,85.85.-|-j,73.63.-b 


Topological currents is one of the well-known physical 
manifestation when a crystalline solid possesses a finite 
Berry curvature 1 The Berry curvature is a geometrical 
property of the Bloch energy band, and acts as an effec¬ 
tive magnetic field in momentum spac^. Hence, topolog¬ 
ical materials may exhibit anomalous Hall-like transverse 
currents in the presence of an applied electric field, in ab¬ 
sence of a magnetic field. In topological insulators®, topo¬ 
logical bands with non-trivial Berry phase leads to propa¬ 
gating surface states that are protected against backscat- 
tering from disorder and impurities. In transition metal 
dichalcogenides, the two valleys carry opposite Berry cur¬ 
vature, or magnetic moment, giving rise to a bulk topo¬ 
logical charge neutral valley current:^ 7 . Recently, these 
bulk topological currents were also experimentally inves¬ 
tigated in other Dirac materials, such as gapped graphene 
and bilayer graphene systenPSl. 

In this letter, we examine the nature of the topological 
currents in black phosphorus (BPJP^HUI when its inver¬ 
sion symmetry is deliberately broken. Unlike more well- 
studied 2D materials such as graphene and transition 
metal dichalcogenides, several key differences are notable. 
First, in monolayer BP, the conduction and valence band 
edges are located at the T point of the rectangular Bril¬ 
louin zone 13 . Second, the energy bands exhibit strong 
anisotropy along its two crystal axes. We will show be¬ 
low that these salient features will lead to a linear trans¬ 
verse neutral topological currents, accompanied also by a 
non-linear transverse charge current. These topological 
currents are maximal (zero) when the in-plane electric 
field is applied along the zigzag (armchair) crystal axes. 
We discuss how these topological effects can be detected 
electrically and optically. 

In this work, we consider monolayer BP, described 
with a four-band tight-binding model 13 as illustrated in 
Fig. la. We construct Bloch-like basis functions |xj) = 
exp[ik ■ (R + by)] |</> R ) where R and bj are the lat¬ 
tice and basis vectors, and the index j runs over all the 
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FIG. 1: (a) Crystal structure of monolayer black phospho¬ 
rus, and the hopping parameters included in the tight-binding 
model. Here, x is along the so-called armchair direction, while 
y is along zigzag, (b) Top view of the crystal structure and 
the primitive unit cell is also indicated, (c) and (d) plots the 
energy dispersion, E k, and its Berry curvature, Hk, for the 
conduction band in the vicinity of the band minimum at T 
point. Inversion symmetry is broken by applying the follow¬ 
ing potentials to the basis atoms, 14,4 = OeV, Vi = —0.5eV, 
and V 3 = 0.5 eV. 

phosphorus atoms within the primitive cell, with a single 
orbital on each atom denoted by | • The Hamiltonian 

matrix can then be constructed from, 

Kiji k) = <* k | n U'i) = Xe ik ' (R+bj ~ bi) /ijj (1) 

R 

where h t j = (<$ | PL |</> R ) contains the tight-binding hop¬ 
ping parameter. The hopping parameters 13 ^ used in this 
work are t\ = —1.22eV, t 2 = — 3.67 eV, £3 = —0.205 eV, 
t -4 = — 0.105 eV, and £5 = —0.055 eU. The secular equa- 
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tion to be solved is Hk | $«k) = E n k | \H„k), where Hk is 
the 4x4 matrix of elements Hij( k), |^„k) and E n k are 
the eigenvectors and eigen-energies. 

The Berry curvature for the electronic Bloch states of 
the n th band can then be computed fronP, 

fj„(k) = fi„(k)z = V k x (^ nk | *V k I’J'nk) (2) 
and it’s magnitude has the following explicit form, 


fi„(k) = 23 


E 

rrij^n 


(^nk| idkx l^mk) (^mk| id ky |\H nk ) 

(Em k - E nk ) 2 


( 3 ) 


Time reversal symmetry implies that Jl„(k) = 
-n„(-k), while crystal lattice with inversion symmetry 
would requires fi n (k) = £l n (— k) = 0. Hence, inversion 
symmetry breaking is necessary to generate a finite Berry 
curvature. 

We consider some basic symmetry properties of the 
Hamiltonian. When the on-site potentials Vj are zero, 
there are two inversion centers i.e. between atom 1 and 
2, and atom 2 and 3. We denote these space inversion 
symmetries as V\ and V 2 respectively. In addition, we 
have mirror symmetry operations A4 X , which interchange 
atoms 1 and 2 with atoms 4 and 3 respectively, and M y 
which is diagonal in the atomic species space. We con¬ 
straint ourselves to inversion symmetry breaking schemes 
with only electrostatic on-site potentials Vj within each 
unit cell. First, consider a perpendicular electric field, i.e. 
Vi = V 2 ^ V 3 = 14 , it breaks M x and V 2 , but not M y 
and V\. Second, consider electrostatic potentials stag¬ 
gering along y (zigzag) direction, i.e. V 2 = V 3 ^ V\ = V4. 
This scheme breaks only Vi, but not V 2 , M. x and M y . 
On the other hand, an electrostatic potential staggering 
along x (armchair) direction will break all inversion sym¬ 
metries, and M x , except for A4 y . Although this configu¬ 
ration can generate non-zero Berry curvature, additional 
symmetries ensure zero transverse currents in this case, 
a subtle point we will elaborate below. 

Alternatively, the combination of an out-of-plane elec¬ 
tric field and an in-plane lattice commensurate electric 
field directed along y direction can break all inversion 
symmetries, since V 4 — V 3 = V\ — V2 0. For ex¬ 
ample, let’s consider V\a = OeV, V 2 = — 0.5 eH, and 
V 3 = 0.5 eV. Fig. lc-d plots the energy dispersion, Ek, 
and its Berry curvature, flk, for the conduction band 
in the vicinity of the band minimum at T point. The 
smaller electron effective mass along the armchair direc¬ 
tion, x, leads to stronger dispersion as shown. Indeed, 
fik is finite, and can take either signs in the Fermi sea 
of the r valley. This is reminescent of the conventional 
valley physics, e.g. in gapped graphene and transition 
metal dichalcogenides where the Berry curvature of the 
two valleys bear opposite signs, except we have only sin¬ 
gle valley in this case. 

In the above-mentioned scheme, the mirror symmetry 
A 4 y is maintained, even when its inversion symmetries 
are being deliberately broken. Time reversal symme¬ 
try requires that its energy dispersion respects E n ( k) = 
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FIG. 2: (a) ‘Bulk Fermi sea’ contribution to the transverse 
current, which can be partitioned into the forward/backward 
propagating currents J^i k - These current are linear with 
electric field £, and persists even in equilibrium, (b) ‘Fermi 
surface’ contribution, also be partitioned into Jtur ; has a £ 2 
dependence, and is a non-equilibrium phenomenon, (c) Il¬ 
lustration of the various transverse current components in a 
typical two terminal device, (d) and (e) plots the dependence 
of the bulk and surface transverse currents with Fermi level 
y, expressed in terms of their conductivities (see text). Cal¬ 
culations assumed temperature of 10 K. 


E n (— k). The above-mentioned mirror symmetry would 
entails E n (k x ,k y ) = E n (k x , —k y ). These symmetries in 
combination also implies E n (k x ,k y ) = E n (—k x ,k y ). For 
the Berry curvature, which is anolagous to an effective 
magnetic field, to produce physical observable that re¬ 
spect this A4 y mirror symmetry would require Vt n (k y ) = 
— Q n (—k y ). Time reversal symmetry would then impose 
the additional constraint that fl n (k x ) = Cl n (—k x ). This 
accounts for the form of the computed Berry curvature 
shown in Fig. Id. 

Following these considerations, we examine the nature 
of topological currents arising from the induced Berry 
curvature, within the semiclassical Boltzmann transport 
theory. In the presence of an external electric field £, the 
carrier velocity acquires a non-classical transverse term 
due to Berry curvature as showrP, 

v„(k) = yVk-EVik - x fi n (k) (4) 

n a 
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The transverse currents can be partitioned into contribu¬ 
tions from the forward and backward propagating states, 
denoted by J ± . By the former (latter), we refer to states 
whose v ra (k) is such that S = sign(v n • £) = ± respec¬ 
tively. J ± can be computed semiclassically up to second 
order in the electric field, i.e. <r^ ;fe £ + af urf l; 2 + 0(£ 3 ), 
with the conductivities defined as, 

a bulk = ~l[ dk/o(k)tt(k) 

°trf = f rfk[V k /o(k) • mV (5) 

where /o(k) is the Fermi-Dirac distribution function, and 
r is the electron scattering time. The linear contribu¬ 
tion to J is a ‘bulk Fermi sea’ phenomenon, and can be 
partitioned into the forward/backward propagating cur¬ 
rents J^ ik as illustrated in Fig. 2a. These are counter- 
propagating currents which persist even when the system 
is in equilibrium. The transverse currents are certainly 
charge neutral since J^ uik = Nevertheless, recent 

experiments have shown that such transport effec ts c an 
be detected via a non-local transport measurementP^. 

On the other hand, the nonlinear contribution to J is a 
‘Fermi surface’ contribution, and can also be partitioned 
into J^ r , see Fig. 2b. These currents, however, are non¬ 
equilibrium in nature, and requires an electrochemical 
potential bias. For example, Fig. 2c illustrates the flows 
of the various current components in a typical two termi¬ 
nal device under bias. Unlike the bulk currents, the latter 
has a net charge current since J+ ur = J^ ur - We contrast 
this with conventional valley physics where Fermi surface 
contributions lead to transverse charge neutral currents 
instead, due to exact cancellation from the two valleys. 
Recent work found that finite nonlinear current can arise 
when the two valleys are not isotropic^, leading to par¬ 
tial cancellation. In a single valley system like BP, the 
effect will be maximal. 

The symmetry of the Berry curvature, i.e. f2„(fc y ) = 
— Cl n (~k y ) and f2„(fc x ) = f2 n (—k x ), has important con¬ 
sequences on the orientation dependence of the various 
transverse currents. J^j k and J^ r will be zero when the 
electric field £ is directed along the armchair direction. 
These transverse currents attained their maximum when 
£ is directed along zigzag direction. 

Fig. 2d-e plots the dependence of these bulk and sur¬ 
face transverse currents with Fermi level /q expressed 
in terms of their conductivities. In 2D electron gas, 
q oc n, where n being the carrier density. We found 
that Ubuik cx fj, and cr sur f oc This is consistent with 
the fact that the number of bulk and surface Fermi states 
scales with q and yjji respectively. The computed linear 
response is an order of magnitude smaller than that pre¬ 
dicted for transition metal dichalcogenideiPl However, 
the finite nonlinear response in this case can produce a 
comparable or larger effect at higher driving fields. 

To gain deeper insights into the above-mentioned is¬ 
sues, we consider an edge-free model system, an electro¬ 
static junction with a built-in electric field aligned along 
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FIG. 3: (a) Schematic illustrating the electrostatic junction 
at equilibrium, along the zigzag y direction. Fermi energy 
q and built-in potential V are indicated, (b) Distribution of 
transverse current density flowing along x , plotted across 
the junction, under different applied out-of-plane potential A, 
and crystal potential A y as indicated (see text), (c) and (d) 
show the integrated transverse current under different applied 
V and Ay respectively. All calculations assumed q = 0.5 eV 
and zero temperature. 


the zigzag (i.e. y ) direction. The electrostatic junction 
provides a built-in electric field that drives the transverse 
bulk current response. Here, we consider the device un¬ 
der equilibrium condition. 

Our approach solves the quantum mechanical scatter¬ 
ing problem microscopically within the above-mentioned 
tight-binding model. We describe this built-in junction 
with V(y) = Utanh(^), as shown in Fig. 3a, where V is 
the built-in potential. In this calculation, we assumed 
a junction transition length of a = 60 nm. The Fermi 
energy, q, is biased within the lowest conduction band of 
BP. We consider the combined effect of an out-of-plane 
electric field, V\ = Vi = ^A z and V 3 = Vi = —^A z , 
and a periodic crystal potential commensurate with the 
BP unit cell, Vj = V 4 = \A V and V 2 = V 3 = —^A y . 
Periodic boundary condition is imposed along x, hence 
k x can be regarded as a good quantum number. Assum¬ 
ing semi-infinite perfect leads, then what we have is es¬ 
sentially a one-dimensional quantum scattering problem. 
The scattering wave functions of the tight-binding prob¬ 
lem, as well as key local observable (e.g. current density, 
the quantity of interest here) , can be solved numerically 
through standard approaches 15 16 . In the following, we 
show that a transverse current can indeed flows along the 
junction, i.e. x direction, when appropriate symmetries 
are broken in consistent with the semiclassical discussions 
above. 

We are interested in the transverse currents arising 
from the forward propagating states deep in the Fermi 
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FIG. 4: Schematic illustrating the excitations of states with 
positive or negative Berry curvatures with circularly polarized 
light. For example, right circularly polarized light couples to 
the states with positive Berry curvatures, which will produce 
an electrical current flowing towards the right contact and the 
top edge. The transversal currents changes sign with A z . 


sea, as depicted by the shaded part of the energy disper¬ 
sion in Fig. 3a. Here, all calculations assumed [i = 0.5 eV, 
and zero temperature. In Fig. 3b, we plot the calculated 
transverse current density profile ( y ) across the junc¬ 
tion, assuming V = 0.2 eV. When A y and A- are non¬ 
zero, crystalline inversion symmetry is broken and a finite 
Berry curvature is present. Indeed, a finite transverse 
current can be observed, which peaks at the middle of 
the junction where electric field is maximum. Away from 
the junction, goes to zero. 

Symmetry arguments inform us that the transverse 
current should be odd with respect to V i.e. V) = 

—JxW)- The Hamiltonian describing the device with 
the built-in junction is related by the mirror symme¬ 
try M y , in the small V limit. Since M y does not af¬ 
fect the transverse current, we have J x [—V) = J.J (V). 
On the other hand, time reversal symmetry would re¬ 
quire J X (V) = — Jx(y)- Taken together, we have 
J x {—V) = —J X (V). Hence, the response to a small 
applied V should be linear. Numerical results shown in 
Fig. 3c confirms this simple argument. 

From the numerics, we found that the transverse cur¬ 
rent have two distinct contributions, i.e. Z+ ~ AJ + 


A y A z £. These trends can be observed in Fig. 3c and d. 
The latter contribution is analogous to the ‘bulk Fermi 
sea’ semiclassical Jbuik we discussed earlier 17 . At small 
V, we clearly distinguish the linear regime consistent 
with the semiclassical result Jbuik oc £. We observed a 
rollover in Z+ at larger positive V, as more states deep in 
the Fermi sea are completely reflected due to the energy 
barrier. When V is negative, there is no energy barrier, 
hence the linear trend persists. The former contribu¬ 
tion to X x is a residual transverse current that probably 
has a different origin. A finite A 2 breaks M x symmetry 
and produces a current along x. Since the magnitude 
of Z+ is tunable by varying A z , it provides an obvious 
way for the detection of the proposed effect. The pro¬ 
posed effect can also be detected optoelectronically. Fig. 
4 illustrates a possible experimental scheme. Circularly 
polarized light can couples preferentially to states with 
Berry curvatures of particular sign, and producing a lon¬ 
gitudinal and transverse electrical current. 

Last but not least, we discuss some considerations on 
the experimental observation of this effect. Breaking of 
the crystal inversion symmetry is key. Finding the appro¬ 
priate substrate which is commensurate along the zigzag 
direction of BP is needed to provide a finite A y . Cer¬ 
tainly, high mobility samples are desirable for the obser¬ 
vation of the proposed effect. Recently, encapsulation 
of BP with hexagonal boron nitride 1 2 3 4 5 6 -^, all within a con¬ 
trolled inert atmosphere, has allowed for higher carrier 

mobilitie^HUl 

Indeed, high quality BP has made possi¬ 
ble the first observation of prominent quantum magneto¬ 
oscillations in these devices!®®!. The results we obtained 
here are also applicable to other emerging 2D materials 
with broken i nver sion symmetry and anisotropic bands 
at the r valley®®^®. 
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